34th 



International Conference on High Energy Physics, Philadelphia, 2008 




mixing at DO experiment 



T. Moulik (for the DO collaboration) 

Department of Physics, University of Kansas, Lawrence, KS 66045, USA 



In this report, we present a preliminary measurement of the mixing parameter using samples of four partially 
reconstructed semileptonic decays and one fully reconstructed hadronic decay mode, corresponding to approx- 
imately 2.4 fb~^ of integrated luminosity. We perform an unbinned likelihood fit to the proper decay length and 
obtain Arris = 18.53 ± 0.93(stat.) ±0.30(syst.) ps'^. 



1. INTRODUCTION 



Particle-antiparticle oscillations are observed and well established in the system. The mass difference Arud is 
measured to be Arrid = 0.507 ± 0.005 ps""'^ mesons are known to oscillate with a high frequency according 

to standard model predictions. Observing the oscillations in the system has been an important focus of the 
B physics program at both the DO and the CDF experiments at Tevatron. DO reported direct limits on the B^ 
mixing parameter ArUg Q using the B'^ D7 fJ.'^vX, D~ — > (j)w~ decay mode ^. CDF reported a measurement 
of this parameter exceeding 5 a significance y. The measurement of this parameter is an important test of the 
CKM (Cabibbo Kobayashi Maskawa) formalism of the standard model, and combining it with a measurement of 
Anid allows us to reduce the error on Vtd and constrain one side of the CKM triangle. This report describes the 
measurement of the 5° mixing parameter at the DO experiement using 2.4 fb~^ of data with 3 a significance. 

We use the central tracker, muon chambers and calorimeters to reconstruct the B decays. Details of the detector 
can be found elsewhere [3]. We use a single inclusive muon trigger or a di-muon trigger to accumulate samples for 
mixing studies. In the case of di-muon trigger the other muon acts as the tag muon used to identify the flavor 
of the B meson which we discuss later. The trigger requires a good muon identified by the muon chamber with a 
matching track in the central tracker in the pseudo-rapidity range of |?7| < 2.0. The triggers use pt cuts between 
3 — 5 GeV/c and the trigger is prescaled or turned off depending on the luminosity. Hence for the i?° — * D'^e^vX 
decay mode, and the hadronic decay mode, we are using a tagged sample with the muon acting as a tag. 



DO reported direct limits on B^ oscillations using the decay mode i?" D^fi^iyX decays with 07r~ 2J. 

In this report, we present a preliminary measurement of the -B" oscillations [5], using four additional decay modes. 
We use three additional semileptonic decay modes namely Dje^vX, with Dj (pir' and i?° Dj fi^vX 

decays with D~ — > K*'^K~ , K*^ if+7r~ (The if*" and the (j) candidates are required to be consistent with known 
mass and width [l| of these two resonances) and KgK^ . We also use a fully reconstructed hadronic mode, 

namely, B^ D'^tt^ + X decay, with (pn^ . 

B1 decays are identified in their semileptonic modes, both in muon and electron mode. Muons are required to 
have Pt > 1.5GeV/c, and to be in the pseudo-rapidity region of |r/| < 2.0. Electrons are required to have pt > 2.0 
GeV/c and are identified in a pseudo-rapidity region oi\-q\ < 1. The Dj and B^ decay products are constrained to 
originate from a common vertex and the B^ and D~ decay vertices are required to be significantly displaced from 
the pp collision vertex. 



2. DECAYS SAMPLE SELECTION AND RECONSTRUCTION 



Charge conjugated states are implied throughout the text 
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The K**^ — > K^TT^ decay mode requires special treatment on account of large reflections, as both real physics 
processes and combinatorial background contribute to the signal peak. The peak is comprised of the signal mode, 
L»7 K*"K-, K+n-, the physics processes, D+ K-tt+tt+ or D+ K*^tt+{K*° if+Tr"), A+ 

K^TT^p^, K*^K^{K*^ if+TT^) (Cabibbo suppressed) and combinatorial background. We fit for these 

contributions in an unbinned likelihood fit. Fig. [1] shows the fit to the mass distribution of the AIkttK system with 
the individual contributions superimposed. 




Figure 1: The M(^^)a' invariant mass distribution where Ds —> K K for 5° — > Ds fi^vX decay mode 



3. FLAVOR TAGGING 

The flavor of the initial state of the is determined using a likelihood ratio method, based on the properties of the 
other b-hadron in the event (opposite side tagging) and the properties of the particles produced in assosciation with 
the reconstructed i?^ meson (same side tagging) and we then combine the OST and the SST taggers. The performance 
of any tagger is characterized by its efficiency defined as £ = N tag /N tot (where Ntag is the number of tagged 
mesons and Ntot is the total reconstructed B° mesons), and the dilution V which is defined as I? = Nrs^^Nws (^-^j^ere 
Njis and N^ys are the right-sign and wrong-sign tags respectively). 

The OST is calibrated using B —^ ^i+uD*^ data events, and we obtain the dilution 23 as a function of a tag variable d 
(whose sign indicates a 6 or 6, and its value indicates the "6"-ness of the tag), to provide an event-by-event "predicted" 
dilution which is used in the unbinned likelihood fit described in section HI More details on the development of the 
OST and results can be found in @. The total OST effective tagging power of eV^ = (2.48 ± 0.2i;° °^)%. 

The combined tagger is calibrated using decay modes, — > J/tp(p (using both data and MC) and 5° ^Ds{(f)'K) 
(using MC only) (See Fig. [1). The total combined tagging power is eV^ = 4.49 ± 0.88%. More details on the SST 
and the combined tagger can be found in [7] . 

The SST is used by i?g pi^vX, with (j>TT~ decay mode only. All the other decay modes use the OST. 



4. UNBINNED LIKELIHOOD FIT AND AMPLITUDE SCAN 
4.1. Visible Proper decay length and scale factor 

In semileptonic decays the proper time gets smeared due to the presence of neutrino. To take this into account we 
introduce a K factor estimated from Monte Carlo (MC) simulation. It is defined as K — pt{IDs)/pt{B). The real 
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Figure 2: The Combined Tagger dilution calibration 



proper decay length is related to the measured or visible proper decay length (Im), by the relation ct{B'^) — ImK 
where Im — M[B^) ■ (Lt)/ {pt{IDs)) is the measured visible proper decay length (Im)- Lt is the distance from the 
primary vertex to the decay vertex in the transverse plane projected onto the IDs momentum, and M(B^) is the 
mass of the 5° meson as obtained from PDG ^l^. 

The Im uncertainty is determined by the vertex fitting procedure, track parameters, and track parameter un- 
certainties. To account for any imperfections in modeling of detector uncertainties, we use scale factors that are 
dependent on the topology of hits in the silicon micro-strip tracker to correct the Im uncertainty. These scale factor 
corrections were obtained from QCD data and MC samples. 



4.2. Likelihood Fit 

An unbinned likelihood fit is used to describe and fit for the oscillation. All flavor tagged events with 1.72 < 
M(Ds) < 2.22 GeV/c^ are used in the fit. The probability for oscillated (osc) and non-oscillated (nos) B^ decays, 
as a function of true proper decay length, x, can be writtten as 

= — • e"^ • 0.5 • (1 + 2?cos(Ato, • Kx/c)) (1) 

CTb, 

pT{x) = • • 0.5 ■ (1 - X'cos(Am, • Kx/c)) (2) 

ctb. 

The sample is mostly composed of -B^ decays with some contributions coming from Bd and Bu mesons also. Similar 
equations would hold for these decay modes with the oscillatory term getting modified appropriately. 

Taking into account momentum uncertainty, and convoluting over detector resolution, the probability for the j-th 
decay mode, as a function of the measured proper decay length, can be written as below : 



nos/ osc 



. SWA-)....,,.",/.. „(. .)pr-'-(.,A-) (3) 



fj{K) is the K-factor distribution and we sum over the K— factor bins. ej{x^) is the reconstruction efficiency as 
a function of lifetime cuts. The fraction of each decay mode contribution is calculated from Monte Carlo. A similar 
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probablity equation holds for the combinatorial background {pbg)- The background contribution comes from prompt 
decays with the Djfj,{e) vertex coinciding with the primary vertex, background with fake vertices distributed around 
the primary vertex, and long lived background. 
The likelihood can then be written as follows : 

iV 

= - fsig)Pi,bg + fsigPi,stg) (4) 

i=l 

The amplitude method was first proposed elsewhere [s!]. It involves modifying the likelihood term by introducing 
an amplitude term A in front of the oscillatory cosine term, such that. 



Cocl±AT> cos{Amst) 



(5) 



We minimize —2lnC fitting for the parameter A while T) is known and Am^ is varied. The value of Am^ where A 
is consistent with 1 and inconsistent with then gives the measurement of the Am^ parameter. All values of Atos 
for which A+ 1.645 < 1 are excluded at 95% confidence level. The sensitivity of the mixing measurement is 
defined as the Amg value for which 1.645 aj^ — l. 

The combined amplitude scan and the difference between log-likelihood at minimum and elsewhere (Aln£), as a 
frmction of Am^, can be seen in Fig. [3] Aln/^ is derived from the amplitude scan. More details can be found in 8]. 
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Figure 3: The combined amplitude and likelihood scan. 



From a parabolic fit around the minimum of the Aln£ distribution, we obtain Amg = 18.53 with statistical 
error of 0.93 ps^^. The systematic error is estimated as 0.3 ps^^ , thus giving us a preliminary measurement of 
Arris = 18.53±0.98ps^^ (stat. -fsyst.). From the depth of the likelihood minima, we calculate an overall significance 
of this measurement as 2.9 a. 
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